Bias- and bath-mediated pairing of particles driven through a quiescent medium 
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A particle driven by an external force in a molecular crowding environment - a quiescent bath of 
other particles, makes their spatial distribution inhomogeneous: the bath particles accumulate in 
front of the biased particle (BP) and are depleted behind. In fact, a BP travels together with the 
inhomogeneity it creates. A natural question is what will happen with two BPs when they appear 
sufficiently close to each other such that the inhomogeneities around each of them start to interfere? 
In quest for the answer we examine here, via Monte Carlo simulations, the dynamics of two BPs in 
a lattice gas of bath particles. We observe that for a sufficiently dense medium, surprisingly, both 
BPs spend most of the time together which signifies that the interference of the microstructural 
inhomogeneities results in effectively attractive interactions between them. Such statistical pairing 
of BPs minimizes the size of the inhomogeneity and hence reduces the frictional drag force exerted 
on the BPs by the medium. As a result, in some configurations the center-of-mass of a pair of 
BPs propagates faster than a single isolated BP. These jamming-induced forces are very different 
from fundamental physical interactions, exist only in presence of an external force, and require the 
presence of a quiescent bath to mediate the interactions between the driven particles. 



I. INTRODUCTION 

A biased particle (BP) traveling in a bath of parti- 
cles, which move randomly without any preferential di- 
rection, drives their spatial distribution out of equilib- 
rium. The bath particles accumulate, creating a "traffic 
jam" in front of the BP and are depleted behind it. This 
BP can be, e.g., a charge carrier subject to an electric 
field or a colloid moved with an optical tweezer. The 
bath particle may be, e.g., colloids dispersed in a solvent 
or adatoms performing activated hopping motion among 
the adsorption sites on a solid surface. 

Such microstructural changes, which substantially en- 
hance the drag force exerted on the BP, have been ob- 
served experimentally; in particular, in microrheological 
measurements of the drag force on a single colloid driven 
through a A-DNA solution [l| or for a biased motion of 
an intruder dragged into a monolayer of vibrated grains 
0. Formation of an inhomogeneous nonequilibrium dis- 
tribution has also been revealed by Brownian Dynamics 
simulations of a driven colloid in a A-DNA solution [I], Q 
and in a colloidal crystal In the latter case, it was 
shown that a large enough BP generates a sufficient stress 
to produce defects, which remain localized near the BP 
and affect the frictional drag force. 

Microstructural changes of a quiescent medium caused 
by a biased probe were extensively studied analytically 
[oT-[ll| for hard-core lattice gases with simple exclusion 



dynamics, in which all particles except one have sym- 
metric hopping probabilities, while one of them - the BP 
- has a preferential direction of motion. 

In one-dimensional systems the size of the jammed re- 
gion in front of a BP (as well as the size of the depleted 
region in the wake) grows in proportion to the traveled 
distance. Thus the jamming- induced contribution to the 
frictional drag force 7 exerted on the BP by the bath 
particles exhibits an unbounded growth, 7 ~ t 1 ! 2 , (t be- 
ing time), so that the BP velocity V t vanishes]!, [fjj|, 
V t (1) ~ t- 1 ' 2 as t ->■ 00. This insures the validity 
of the Einstein relation for anomalous diffusion in onc- 
dimensional hard-core lattice gases d, 0] ■ 

In higher dimensions, the BP velocity attains a drift 
value = VfrZ^ and the bath particle distribution 

reaches a non-equilibrium stationary form [sMTlj . The 
density profiles are strongly anisotropic with a traffic- 
jam like region in front of the BP and a depleted region 
in its wake. Strikingly, behind the BP the bath particle 
density approaches the mean value p as a power-law of 
the distance x: \jx z l 2 and ln(a;)/izi 2 in 2D and 3D [1- 
[Tlj . which signifies that the medium "remembers" the 
passage of the BP on large temporal and spatial scales. 
The drift velocity V^ 1 ) and the jamming- induced drag 
force 7 have been determined for the lattice gas model 
7f 
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and also for a driven probe in a colloidal mixture 
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The next step in the understandin g of the jamming- 
induced forces has been done in refs. [14 - llq |. Dzubiella 
et al [HI studied the effective interactions between two 
fixed colloids in a quiescent viscous solvent exposed to a 
flowing bath of small Brownian particles, while Krueger 
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and Rauscher and Khair and Brady [l6| considered 
the case of two colloids translating along their lines of 
centres with fixed velocities and a fixed distance apart of 
each other in an otherwise quiescent dispersion of nonin- 
teracting colloids. It was realized that, remarkably, mi- 
crostructural changes induce effective forces between two 
colloids, which may be either repulsive or attractive, de- 
pending on their mutual orientation. 

In this paper we pose a very natural question within 
the context of microfluidics/microrheology or biased dy- 
namics under molecular crowding conditions: What will 
happen with two BPs in a quiescent medium of mutu- 
ally interacting particles when both, move not with a 
prescribed velocity along some fixed lines, but rather 
perform a biased random motion subject to some ex- 
ternal force? In contrast to the situations studied in 
refs. [1, [3 [HI , here the BPs can change their relative 
position in space and hence, by monitoring their trajec- 
tories we can understand the overall effect of microstruc- 
tural changes of the medium on the interactions between 
them. 

In order to be as transparent as possible, here we resort 
to a minimal model of a hard-core lattice gas of parti- 
cles whose dynamics obeys the so-called simple exclusion 
process (SEP) [13, dl- We note parenthetically that this 
model of dynamics is quite realistic and applies to many 
physical systems, such as.ejr dynamics of adatoms on 
solid surfaces (see Refs. |8l-[li1. [iH HH for more details 
and other systems). We note, as well, that it allows us to 
single out the effect of microstructural changes, i.e., bath- 
mediated interactions, and to exclude possible effects of 
solvent (if present) and solvent-mediated interactions [l9| 
between the BPs. 

In this model the lattice gas particles - the bath parti- 
cles - have symmetric hopping probabilities while two 
particles - the BPs - are subject to an external force 
and have asymmetric hopping probabilities. Tracking the 
BPs trajectories in Monte Carlo simulations, we observe 
a phenomenon of statistical pairing of biased particles. 
We realize that for sufficiently dense systems the frac- 
tion of time which the second BP spends at a given point 
in space has an apparent maximum in the vicinity of the 
first BP. Hence, the interference of non-equilibrium den- 
sity profiles of the bath particles formed around each of 
the BPs results in an effective attractive interaction be- 
tween them. Apart of this, we analyze the properties 
of the jamming-induced frictional forces and determine 
the velocities of the BPs appearing in different configu- 
rations. 

The paper is outlined as follows: in Section 2 we de- 
fine the model. In Section 3 we focus on dynamics of a 
single biased particle in a quiescent bath, describe the 
density profiles of the bath particles around the BP and 
determine the force-velocity relation. In Section 4 we 
consider dynamics of two BPs. Here we define most prob- 
able paths that the BPs follow and describe the density 
profiles of the bath particles forming around a pair of 
BPs appearing in different configurations. Apart of this, 




FIG. 1. A lattice with randomly moving hard-core 
particles. Bath (blue) particles have symmetric hopping 
probabilities. BPs (red) are subject to a constant force and 
have asymmetric hopping probabilities. 



we define the frictional force exerted by the medium on 
the BPs and determine the velocity of a pair of BPs. We 
conclude in Section 5 with a brief recapitulation of our 
results and an outlook of future work. 



II. THE MODEL 

Consider a square lattice of S = L x x L y sites, of 
spacing a and with periodic boundary conditions. The 
lattice is populated with two different types of parti- 
cles: N - M bath particles and M (M = 1,2) BPs 
(see figure [T]). All particles experience hard-core inter- 
actions, such that each site can be cither empty or oc- 
cupied by at most one particle. The state of each site 
(X, Y) is described by a time-dependent occupation vari- 
able T)(X,Y); T]{X,Y) = 1 if the site (X, Y) is occupied 
and n(X, Y) — 0, otherwise. 

Particle dynamics is defined by the so-called SEP - 
simple exclusion process [13, EU • Each particle bears an 
exponential clock; in general, the mean jump-time of the 
bath particles and that of the BPs may be different; say, 
it is t* for the bath particles, while for the BPs it is 
t. We will focus in what follows on the simplest case 
t* = t = 1 and will only briefly mention the effects of 
different jump times on the friction coefficient. 

When the clock rings, a particle attempts to jump from 
the site (X, Y) it occupies to one of the four nearest- 
neighboring sites (X' , Y') according to the normalized set 
of probabilities p{X, Y\X', Y'). Once a jump direction is 
chosen, and if the destination site is empty, the particle 
moves to it, otherwise it remains on the site it occupies. 
This stochastic exclusion dynamics is a Markov process 
on a state space M = {0, 1} L * ® {0, 1} L «. 



3 




We stipulate next that the dynamics of the bath par- 
ticles obeys a symmetric SEP (all hopping probabilities 
= 1/4), while the BPs are driven by an external field 
F = —Fix and evolve according to an asymmetric SEP: 

p(r, r ± ae x ) = Z~ 1 e^ (TF/2 , p(r, r ± ae Y ) = Z' 1 , 

(1) 

where ex — (1, 0) and ey = (0, 1) are unit shift vectors, 
r = (X,Y), Z = 2(1 + cosh (^crF/2)) and f) is the inverse 
temperature. 

III. ONE BIASED PARTICLE 



To set up the scene, we focus first on the case of a 
single BP on a lattice with N — 1 bath particles. As we 
have already remarked, a single BP produces microstruc- 
tural changes in the medium it travels in, driving the 
spatial distribution of the bath particles out of equilib- 
rium. To quantify the microstructural changes of the 
medium, we consider the following realization-dependent 
" inhomogencity" measure: 



9o (r;t) = 



\J*dT^J2S(T i (r)-B. 1 (T)-T)\ , (2) 



where S(-) is the Kronecker-delta, which equals 1 when 
its argument is and is zero otherwise, p = (N — 1)/(S — 
1) is the mean density of bath particles, while Ri(t) 
and rj(r) denote the positions of the BP and of the i- 
th bath particle at time r, respectively, for a given set 
of realizations of their trajectories. An average of go(r; t) 
over different realizations of particles' trajectories defines 
the time-averaged van Hove function |20[. 

The realization-dependent functional go(r;t) in Eqf2] 
defines the fraction of time the site r = (x,y), in the 
frame of reference moving with the BP, is occupied by 
bath particles during a time t for a given realization of 
trajectories of the bath particles and the BP. [? ] If the 
spatial distribution of the bath particles converges to a 
stationary form, i.e., if limt_>oo <7o( r ) — 9o( r ) exists, 
then pgo(r) can also be thought of as the bath parti- 
cles density profile as seen from stationary moving BP. 
Clearly, go(0) = due to the hard-core exclusion and 
go(r) — > 1 when |r| — > oo. We have conveniently normal- 
ized (?o( r ; t) to p, so that any deviation go ^ 1 indicates a 
non zero dynamical correlation between the BP and the 
medium. 

In figure [5] we depict the bath particles density go (r) 
using a colour map. The density profiles around a sta- 
tionary moving BP are characterized by a jammed, high- 
density region in front of the BP and a pronounced region 
depleted by the bath particles past the BP. This agrees 
quite well with the theoretical prediction of Refs.|7H10|. 
We verify, as well, the theoretical prediction that the den- 
sity past the BP approaches the average value p not ex- 
ponentially with the distance x, but as a slow power-law 
. Moreover, in figure [5] we superimpose the average 



FIG. 2. Microstructural changes of the medium pro- 
duced by a single BP. The profile go(x, y) in Eq[2]is shown 
for a lattice comprising 61 x 21 sites at density p = 1/2 and 
/3aF — 5. The vector field shows the average velocity field 
of the bath particles defined in the reference frame of the BP 
(a black square). Velocity vectors with magnitude less than 
0.0005 are not plotted. 



velocity field, which shows that the BP induces a reg- 
ular global motion of the bath particles predominantly 
towards the regions with lower density. Note that sim- 
ilar motion patterns have been observed experimentally 
for biased motion in granular media Q. 

Next, in figure [3] we depict the dependence of the BP's 
drift velocity V^ 1 ' on the applied force F, i.e. the force- 
velocity relation, for p = 1/2. 

The notable feature of the observed force-velocity 
curve is that in the limit of sufficiently small forcing 
(sufficiently high temperatures) the drift velocity shows 
a Stokesian, linear dependence on the applied force, 
— F/£, which signifies that in this limit the fric- 
tional force exerted by the bath particles on the BP is 
viscous. 

According to Refs.@-Q3 (see also Section VI in 
Ref.[2l|), in this linear regime the friction coefficient £ 
can be expressed as a sum of two contributions, 



6 = £,rnf + £,c 



where the first term, 



4r 



0(7* (1-p)' 



(3) 



(4) 



is essentially a mean-field result corresponding to a per- 
fectly stirred monolayer; one may interpret (1 — p)/r just 
as the frequency of the BPs "successful" jump events. 
The second term is a "jamming-induced" contribution: 



At* 



(tt - 2)p 



/3a 2 (l-p) \ + {1- P )t*/t 



(5) 



-3/2 



stemming out of a cooperative behavior in the monolayer 
- a non-linear interplay between the BP dynamics and 
the formation of non-equilibrium density profiles around 
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FIG. 3. Force-velocity relation. Drift velocity V^ 1 ' of a 
single BP vs F (solid squares). The dashed line is the theoret- 
ical prediction V™' = F/£ in which the friction coefficient £ 
is given by Eq(3] Open circles denote the results of the Monte 
Carlo simulations for the drift velocity of a pair of BPs in the 
(2,0) configuration (see the explanations in Section 4.2). 



it. Analogous results for £ have been obtained for three- 
dimensional [llj and one-dimensional systems; in the 
latter case £ diverges as t — > oo. 



Dividing £ m/ by £ coo „, we have 



^COOP 



1 



P 



1 T 

1>~* 



(6) 



One notices that the jamming-induced contribution to 
the friction coefficient dominates when r* 3> r and the 
bath particles mean density is not too small. Conversely, 
the mean-field contribution is clearly the dominant one 
when p C 1 or t* C t. For moderate densities, £ m / and 
Ccoop are comparable. 

Finally, in figure 3] we compare the result in Eq(3] 
against Monte Carlo simulation results for the slope of 
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FIG. 4. Friction coefficient £ versus density. The dashed 
line is the theoretical prediction in Eq[31 while the symbols 
define the numerical simulations results for <j = f3 = 1 and 
t = t* = 1. 



the force-velocity relation in the limit of small forcing, 
at different densities p and r = r* = 1. One observes 
a very good agreement between a theoretical prediction 
and numerical data. 



IV. TWO BIASED PARTICLES 

We turn now to the case of primary interest - two BPs. 
We suppose that initially the BPs are placed at sites 
(X,Y) and (X',Y r ), and that X < X'; following the 
terminology of Khair and Brady [l6| . we then refer to 
the particle initially at (X, Y) as the " leading" BP, and 
the one at (X', Y') - the "trailing" BP (see figured]). We 
recall that the coordinate system defined in the reference 
frame of the leading BP is be denoted by (x,y). 



A. Mutual orientation of the BPs 

We concentrate first on the analysis of the most proba- 
ble paths of the trailing BP in the reference frame of the 
leading BP. For this purpose, we study numerically the 
behaviour of a realization-dependent functional 



1 



t-to 



drS (R 2 (r) -Ri(r) -r) 



(7) 



to 



where Ri(t) and R2(t) stand for the instantaneous po- 
sitions of the leading and trailing particles, respectively, 
for a given realization of their trajectories. Similarly to 
the functional defined in EqJ^l g(r;t) defines, for a given 
realization of the leading and trailing BPs trajectories, 
the fraction of time during the time interval t — to that 
the site r, (in the frame of reference moving with the 
leading BP), has been occupied by the trailing BP. 

In simulations the leading BP is initially placed at the 
origin, the trailing one is placed at position R 2 (0), while 
the bath particles are placed at random, with mean den- 
sity p, on the rest of the lattice. We let the system 
evolve (the bath particles follow a standard symmetric 
SEP while the BPs dynamics obeys an asymmetric SEP) 
for 10 8 time steps, until the density profile of the bath 
particles around two BPs attains a stationary form. After 
this transient period, we define the moment to when the 
trailing BP re-appears again at site R2(0). Then, during 
the next t — t = 10 6 time steps we evaluate g(r;t) by 
recording the number of times each site r in the frame 
of reference of the leading BP has been visited by the 
trailing BP within this realization of the process. 

In figure [5] we plot the average local occupation times 
<g(r;t)>[? } for f3aF = 5 and R 2 (0) = (3,4). Our re- 
sults show that when the inhomogencitics around each 
BP do not overlap, both BPs travel almost indepen- 
dently. In this case, the profile of the average occupation 
times around the initial position of the trailing BP is al- 
most symmetric (see Fig. 0a) with a small second max- 
imum just after the leading BP. The overlap becomes 
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FIG. 5. Most probable paths. The profile of average local occupation times <g(r; t)>, Eq[7] for the trailing BP commencing 
at R 2 (0) = (3,4) (little black knot), (3aF = 5 and densities p: a) 1/4, b) 1/2, c) 3/4 and d) 0.9. High values of <g(r;t)> 
(yellow) indicate the most probable paths that the trailing BP follows. 



significant for cither sufficiently high density p, larger 
driving force f3aF or naturally, when the leading and the 
trailing BPs are close enough. For progressively higher 
densities of the medium particles (see Figs|5l&-<i) we ob- 
serve a considerable qualitative change in the form of the 
profile <g(r;t)>: it becomes considerably more asym- 
metric and is characterized by an apparent " bridge" con- 
necting the leading and the trailing BPs. Hence, the 
probability of finding the trailing BP in the vicinity of 
the leading one is getting progressively higher. 

Further on, we focus on a single (very long) trajectory 
of a trailing BP (see Fig. [5]). We let the system evolve 
for 10 8 time steps to ensure that the bath particles dis- 
tribution around the two BPs reaches a stationary form. 
Then, we evaluate <?(r; t) by tracking the trajectory of 
the trailing BP in the frame of reference of the leading 
one. We have checked that for sufficiently large times (in 
our simulations t — 10 8 ), g(r;t) converges to a station- 
ary function g(r), which moreover, is independent of the 
initial state. 

Figure |5] shows the local occupation times g(r;t) of a 
single trajectory together with the relative velocity field 



of the trailing BP with respect to the leading one. The 
circulating structure of this field indicates that the pair is 
statistically stable, presenting a stable direction (:r-axis) 
along which the trailing and the leading BPs " pair" , and 
an unstable direction (j/-axis) along which the trailing 
BP moves away from the leader. 

This signifies that when the density profiles of the bath 
particles emerging around each of the BPs start to in- 
terfere, the jamming- induced frictional force exerted on 
the trailing BP is no longer parallel to the x-axis (di- 
rection of the external force F) but is tilted by some 
angle pointing towards the wake of the leading BP. As 
a consequence, the trailing BP experiences an effective 
attraction towards the leading one such that both driven 
BPs statistically " pair" . From the microscopic dynamics 
viewpoint, this interference modifies locally the statistics 
of successful hops of the trailing BP, increasing the like- 
lihood of hops towards the wake of the leading BP. 

To substantiate this claim, we have computed numeri- 
cally an instantaneous jamming- induced nonequilibrium 
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FIG. 6. Local occupation times of the trailing BP. 

Colour map of the occupation time of the trailing BP given 
by g(r)/(S — 1), for a lattice comprising 31 x 21 sites, p — 1/2 
and /3aF = 5. The superimposed vector field superimposed 
shows the average relative velocity field of the trailing BP 
measured in the reference frame of the leading BP. Velocity 
vectors with magnitude less than 0.001 are not plotted. 



force 7(1*; t) that we define as 

?(r;t)= (^(tJ-eO-f/^tJ + e,-))^. (8) 

i={x,y} 

where 7j(r(t)) are the local occupation variables defined 
in Section II. Clearly, j(t; t) is the force that is felt by 
the trailing BP being at site r with respect to the leading 
BP at time moment t. 

On average, 7*(r) will be different from zero only if the 
distribution of the bath particles is inhomogeneous. In 
general, as follows from figure [5J -f(r) should depend on 
the position of the trailing BP with respect to the leading 
one. In figure [7] we show the time averaged y component 
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FIG. 7. Force exerted on the trailing BP. Time aver- 
aged jamming-induced force experienced by the trailing BP 
at position (0,7) (open circles) and (4,2) (solid squares), in 
the reference frame of leading one, as a function of f3aF. The 
error bars indicate the standard deviation. 



of the force 7 that the trailing BP experiences when it 
is at position (0, 7) (circles) and (4, 2) (squares) with re- 
spect to the leader. In the first case, the trailing BP is 
far enough from the leader so that the net force along y is 
numerically zero, irrespective of the strength of the field. 
On the contrary, in the second case, the trailing BP is 
in the stable basin of attraction and 7^ is negative, indi- 
cating that (in agreement with Fig. [S]), the force pushes 
the trailing BP towards the leading one. For small fields 
7 y grows linearly with f3aF and saturates at larger fields. 
Of course, the amplitude of the force 7(1") depends on the 
density of bath particles, external force F and, naturally, 
on the temperature, which controls the rate of the bath 
particles' migration and thus their capability for smooth- 
ing down the inhomogeneities created by the BPs. 



B. The BPs velocities 

The (mean) drift velocity V"W of a single isolated BP 
is totally determined by /3a F, p and the rate of the bath 
particles' migration]^, [l(|. For two BPs, when they ap- 
pear sufficiently close to each other such that the inho- 
mogeneities around each of them start to interfere, their 
drift velocities and the velocity of their center of mass 
will also depend on their mutual orientation. One may 
expect that only at large mutual separations of the BPs 
their velocities are equal to V^ 1 ' . 

To clarify this issue, we have studied the velocities of 
the leading and trailing BPs at a fixed force /3aF = 5 and 
at a fixed bath-particle density p = 1/2, but for different 
mutual orientations. Six different mutual orientations of 
the BPs, together with the corresponding microstructural 
changes of the medium, are presented in figure [8] In 
this figure we also depict the velocity field of the bath 
particles. 

Further on, in figure [5] we plot the instantaneous drift 
velocities of the leading and trailing BPs as the function 
of their mutual orientation. In the left panels, along the 
horizontal line, each "tick" corresponds to the position 
(x, y) of the trailing BP in the reference frame of the 
leading one. Between each pair of vertical dotted lines, 
we place nine points with fixed y and x varying from 1 
to 9, i.e., (1,2/), (2,y) to (9,y). From the left lower panel 
one indeed concludes that at large separations the veloc- 
ities of the BPs are nearly the same, and coincide with 
the velocity (blue dashed line) of a single isolated 
BP. Conversely, when the BPs are close enough, their ve- 
locities may be very different from each other, as well as 
from VW. 

At a fixed y, the velocity of the trailing BP is always a 
non-monotonic function of its x-coordinate: it is always 
minimal for x — 0, grows abruptly with the x-coordinate 
and then decays towards The leading BP velocity, 

at a fixed y-coordinate of the trailing BP, is a growing 
function of the x-coordinate of the latter. Despite the 
fact that the configuration (1,0) is the most probable 
(see figure EJ), the velocity of the center of mass (solid 
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FIG. 8. Microstructural changes produces by two BPs. Colour map of the mean local occupation times of the bath 
particles {go(x,y)}, for six different configurations of two BPs, defined in the reference frame of the leading BP (black square), 
with p = 1/2 and fioF = 5. The arrows define the vector velocity field. Velocities with magnitude less than 0.005 are not 
plotted. 



blue line in figure [5]) of such a pair is not the largest 
one, which is a bit counter intuitive. As a matter of 
fact, this is the consequence of the hard-core interac- 
tion between the leading and trailing BPs that hinder 
the motion of the latter reducing the average velocity 
of the pair. The blocking effect is no longer present in 
the next configuration (2, 0) for which (together with the 
configuration (3, 0)), the velocity of the center-of-mass is 
the largest. This signifies that such a pair of the BPs 
creates the smallest microstructural changes, and in re- 
sponse, encounters the least jamming-induced frictional 
force, hence the least possible dissipation. 

Curiously enough, the velocity of the center-of-mass 
of such a pair is always higher than the velocity of a 
single isolated BP. This can be seen in figure |31 in which 
we compare numerical force- velocity relations for a single 
BP (solid squares) and a pair in the configuration (0, 2) 
(open circles), for different values of the driving force. 
The center-of-mass velocity is also higher than the veloc- 
ity of a single isolated BP for the configurations (3,0), 
(4, 0) to (9, 0), (2, 1) to (9, 1) and (4, 2) to (9, 2). On the 
contrary, the microstructural changes of the medium in- 
duce the largest frictional force on pairs (0, y), for which 
the BPs are in perpendicular orientation with respect to 
the field direction, such as, e.g., configurations (c) and 



(d) . These configurations are indicated in figure by 
vertical arrows. Due to the symmetry of the problem, 
in these configurations the velocities of the leading and 
trailing BPs coincide, and are smaller than the velocity 
of a single isolated BP. 

To better understand such a behaviour we study a 
measure of the size of the inhomogeneity created by the 
BPs. This configuration-dependent measure that we call 
E, equals the sum of the mean occupation times of the 
bath particles, over the lattice sites on which the absolute 
deviation of (go(x,y)) from the mean value p exceeds a 
certain threshold Q: 

E = X>o(z,J/)> , (9) 

5g>g 

where Sg = | (ga(x, y)) — p\. To look at the correlation 
between the size of the inhomogeneity and the jamming- 
induced interaction among different BPs, we denote E„ 
(n = 1,2), as the inhomogeneity measure ©, computed 
for n BPs. In the left top panel of figure [5] we present 
the ratio E2/2E1 for different configurations of the BPs, 
clearly showing that the largest (smallest) velocities of a 
pair of BPs occur for smallest (largest) values of E2/2E1. 
In particular, for the pair configurations whose center- 
of-mass velocity is larger than the single BP, we obtain 
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FIG. 9. Drift velocities of the BPs in different configuration Left lower panel.- Magnitude of the a;-component of the 
velocity of the leading (open squares) and trailing (solid circles) BPs for different mutual orientations, /3aF = 5 and p — 1/2. 
The blue dashed line defines the velocity V' 1 ' of a single isolated BP, while the solid blue line defines the center-of-mass 
velocity of a pair of BPs. The arrows indicate some of the configurations depicted in figure [S] Left upper panel. - The ratio of 
"effective" inhomogeneities created by two BPs and a single BP (see the text for the explanations). Right panel.- Drift velocity 
of the transversal configurations as a function of the separation y. The inset shows that the drift velocity converges to V' 1 ' 
exponentially. 



E2/2S1 < 1, indicating that the overall size of the in- 
homogeneity produced by the pair of BPs in such con- 
figurations is smaller than the corresponding size of two 
isolated BPs. 

Finally, in the right panel of figure [9] we show that at 
large separations along the y-axis the velocity of the pair 
approaches the velocity of a single isolated BP exponen- 
tially, as is the case of configuration (f), corresponding 
to the orientation (0, 15). 

V. CONCLUSIONS 

To recap, we have addressed here a problem of effec- 
tive, non-equilibrium interactions that emerge between 
two driven probes in a medium of randomly moving 
hard-core (but otherwise noninteracting) particles - a 
quiescent bath. We have shown that for a sufficiently 
dense medium the probes experience an attractive force 
towards each other, leading to their statistical pairing. 
Such a pairing sets in when the non-equilibrium inhomo- 
geneities in the distribution of the bath particles, created 
by the probes, start to interfere. The inhomogeneities 
around each driven probe decay exponentially with the 
distance from the probe, except for the wake of the probe 
in which the decay is algebraic. Consequently, these 
non-equilibrium interactions are anisotropic and typi- 
cally short-ranged, except for the situations when the 
second driven probe appears in the wake of the first one. 

The formation of pairs reduces the overall size of the 
inhomogeneity, minimizing the frictional drag force the 
medium exerts on each probe. As a result, in some con- 



figurations the center-of-mass of a pair propagates faster 
than a single isolated BP. The jamming-induced forces, 
which emerge in the situation under study, are very dif- 
ferent from fundamental physical interactions, exist only 
in presence of an external force, and require the presence 
of a quiescent bath to mediate the interactions between 
the driven intruders. 

We note that our results have been obtained for a 
somewhat simplified model of a non-interacting lattice 
gas with simple exclusion dynamics, which allowed us to 
single out the effect of the jamming- induced interactions. 
This model can be generalized in several directions. First 
of all, one may consider a situation appropriate to a col- 
loidal solution, when some solvent is present. Solvent 
itself will produce long-range hydrodynamic interactions 
between the driven probes and correlate dynamics of the 
bath particles (see Ref.[lj| and references therein). One 
may expect that the pairing effect will become more pro- 
nounced in this case. Second, we have considered the 
case of just two driven probes. It might be interesting to 
study the specific features of pairing in the situation when 
there are many of them - the effect observed here seems 
very much like "an elementary act" for the phenomenon 
of lane formation in partially driven colloids [22|, |23[. Fi- 
nally, we note that we have focused here solely on the case 
when the biased motion of the BPs results from the pres- 
ence of an external force acting on them. For biased mo- 
tion in intracellular media or under molecular crowding 
conditions, it might also be interesting to consider other 
types of biased motion, e.g., the cases of self-propelled 
particles or swimmers. These important situations merit 
further investigation. 
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